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Abstract: We present a new method to search for heavy nuclei sources, on top of background, in the Ultra-High Energy 
Cosmic Ray data [1]. We apply it to the 69 events with energies E > 55 EeV published by the Pierre Auger Collaboration. 
We find a set of events for which the method reconstructs the source near the Virgo galaxy cluster. The probability to 
have a comparable set of events in some background is ~ 0.7%. The reconstructed source is located at ~ 8.5 degrees 
from the active galaxy M87. The probability to reconstruct the source at less than 10 degrees from M87 for data already 
containing a comparable set of events is ~ 0.4%. This may be a hint at the Virgo galaxy cluster as an ultra-high energy 
heavy nuclei source [2, 1]. We discuss the capability of current and near future experiments to test this possibility. Such 
a scenario gives a self-consistent description of the Auger anisotropy and composition data at the highest energies. 
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1 Introduction 

The Pierre Auger Observatory measurements of the Ultra- 
High Energy Cosmic Ray (UHECR) composition are 
compatible with a shift towards heavy nuclei, above ^ 
10^^ eV [3]. The analysis of the muon data from the 
Yakutsk HAS Array also hint at a significant fraction of 
heavy nuclei at the highest energies [4]. On the contrary, 
both HiRes measurements [5] and preliminary results of 
Telescope Array [6] are still consistent with a proton com- 
position. 

Both for protons and iron nuclei, UHECR sources must be 
located in the local Universe, within r < 150Mpc for en- 
ergies E > 6 X 10^^ eV. They should lie within the Large 
Scale Structure (LSS) of galaxy distribution. 

Until now, methods to look for sources have been proposed 
for proton or light nuclei primaries [7, 8]. We present here 
a new method to search for UHE heavy nuclei sources. 
In most regions of the sky, one cannot detect heavy nu- 
clei sources without a better knowledge of the Galactic 
Magnetic Field (GMF) than currently available [9, 10, 1 1]. 
Nonetheless, we show that in some favourable cases, one 
may construct an algorithm to detect some of such sources 
with the present knowledge of the GMF. 

We apply this method to the 69 Auger events published in 
Ref. [12]. We detect a set of events for which the associ- 
ated reconstructed source lies near Virgo. The Virgo clus- 
ter is in principle a good candidate for containing one or 
several UHECR source(s). The probability to have a com- 
parable set of events in some background and reconstruct 



the source in any direction of the sky is ^ 7 x 10^'^. If 
such a set of events already exists in the data, the proba- 
bility to reconstruct the source at less than 10° from M87 
is ^ 4 X 10^'^. The possibility that the detected feature is 
the heavy nuclei image of Virgo would be compatible with 
both the Auger composition and anisotropy at the highest 
energies. However, the present statistics do not allow us to 
conclude firmly. We discuss the ability of current and near 
future experiments to confirm or rule out this possibility. 

We present the method in Section 2. In Section 3, we anal- 
yse the Auger data of Ref. [12] and report the noteworthy 
feature. Section 4 presents a discussion of the results. 

2 IMethod 

For UHE proton sources, events are expected to be de- 
flected on the sky from the sources as l/i?[7,8]. For heavy 
nuclei primaries, images of sources located in most regions 
of the sky would have more complicated shapes, even at 
the highest energies [13, 14, 15, 9, 10, 11]. We find that in 
some GMF models and some regions of the sky, sources of 
nuclei of charge Z can still have at least one of their images 
above ^ 55 EeV which roughly looks like a proton source 
image enlarged by a factor ^ Z [1]. Such images can be 
detected without a better knowledge of the GMF than cur- 
rently available. The following method is optimized to look 
for such favourable heavy nuclei source images in the data. 

Fig. 1 presents a sketch of the method. It starts select- 
ing an event with energy Ei > lO'^^eV, denoted by "1". 
Then, an assumption on the UHECR deflection power V of 
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Figure 1 : Sketch of the method proposed to detect extended 
UHE heavy nuclei sources S (black disk) with "enlarged 
proton-like" images. The source events are represented 
by decreasing energy order, by the magenta filled and red 
open boxes, green open circles and blue crosses. The sector 
(opening angle O and extension R) is highlighted in grey. 
X' axis in red and S' for the reconstructed source. 
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Figure 2: Portion of the sky in Galactic coordinates, show- 
ing the most significant feature detected with our method in 
the Auger data of Ref. [12]. The 13 events in the sector are 
surrounded with magenta circles. Symbols for the Auger 
cosmic rays according to their energies -see text. Black 
disk for the Virgo cluster, blue line for the X' axis and 
thick red cross for the position of the reconstructed source. 



the GMF [8] is made, and only events at angular distances 

< R — 2?/(55EeV) — V/Ei from 1 are considered here- 
after Events with energies £'2, located at angular distances 

< V/ E2 — "D/Ei from 1, are tested by decreasing energy 
order Let us take one of them, denoted "2" in Fig. 1 . 1 and 
2 define the central line of a sector shaped region, high- 
lighted in grey in Fig. 1. The sector has an angular exten- 
sion R and an opening angle 8. Its vertex corresponds to 
1 . Let us define X' as the angular distance to L If the 
number of events in the sector and their correlation coef- 
ficient Con:{X' ,1/E) are respectively larger or equal to N 
and Cmin, there is detection. Otherwise, the next event 2 
in the ordered list is tested, until there is either detection or 
all events 2 have been tested without any detection. In case 
of detection, the source is reconstructed along the axis de- 
fined by 1 and the center of mass of all points in the sector 
(red axis in Fig. 1), by fitting 1/E versus X' as in Ref. [8]. 
This method has four free parameters: Af, Cmin, T> and <d. 
V and <d respectively depend on the regular and turbulent 
GMF contributions to UHECR deflections. 



3 Application to the Auger data 

We apply the method presented in Section 2 to the 69 Auger 
events, with E > 55EeV, published in Ref. [12]. We scan 
them over discretized sets of values for Af, Cmin, 2? and 
e. A/" e {4, 5,..., 69}, C,„i„ e {-1,-0.9,..., 0.9}. Ac- 
cording to the relative contributions to UHECR deflections 
of the regular and turbulent GMF components found by 
Ref. [16], e G {10°, 20°, ..., 80°} should suffice. Ref. [17] 
notes that deflection angles on the sky for 10^" eV pro- 
tons are expected to be ~ 1 — 2.5°. Therefore, as- 
suming iron nuclei, we take V e {26°, 39°, 52°, 65°} x 
lO^'' eV. We confront the data with some random back- 
ground, made of 69 events distributed on the sky accord- 
ing to the Auger exposure. Their energies are distributed 



within four bins: [SSEeYlQi^-^eV], [10l9■^10l'5■9] eV, 
[10i9^1020]eV and above lO^OeV, according to the 69 
events energies. No noteworthy dependence on the distri- 
bution of energies within the bins was found. For the last 
bin, we take a E~'^-^ spectrum, following [18], and a max- 
imum energy E'max = 10^"'^ eV, because the attenuation 
length of heavy nuclei rapidly becomes shorter than a few 
Mpc above this energy [19]. We checked that assuming 
lower values for -Emax only increases the significance of 
the signal detected below. 

An interesting signal is found in the Auger data around 
the event with energy E ~ 142 EeV, which is located at 
~ 34° from the center of Virgo. This event plays the role 
of "1" in the method presented in Fig. 1. For the tested 
ranges of parameters, the lowest probability to reproduce 
the Auger data with the background is found for Af — 13, 
Cmin = 0.6, 2? = 39° X 10^0 eV and 6 == 40°. In the 
data, Con:{X',l/E) ~ 0.66. These 13 detected events are 
plotted in Galactic coordinates and surrounded with ma- 
genta circles in Fig. 2. The filled magenta boxes, red open 
boxes, green open circles and blue crosses respectively cor- 
respond to the arrival directions on the sky of Auger events 
with energies E > 10^° eV, lO^^-^ eY < E < lO^" eV, 
XQiS-S eW<E< lO^^-^ eV and E < lO^^-* eV. 

Since one cannot know a priori the best values for the four 
free parameters of the method, one has to penalize over 
all of them [1]. After penalization, the probability to find 
an at least as significant feature in the background as the 
feature detected in Fig. 2 is ~ 6.6 x 10"'^. The source 
associated with these 13 events (red cross in Fig. 2) is re- 
constructed at ~ 8.5° from M87, near the Virgo cluster 
position (black disk). Even in case Virgo would only con- 
tain one UHE heavy nuclei source, it would shine as an ex- 
tended source on the sky because of magnetic fields inside 
the cluster [20]. Moreover, due to the poor 1/E ordering 
of events expected for heavy nuclei sources, the accuracy 
on the position of the reconstructed source would not be 
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better than ±10°. Therefore, the detected events are com- 
patible with a common emission from Virgo. (Most of) 
the events in the "Cen A region" (-60° < / < -30° and 
0° < 6 < 30°) and the 142 EeV event may be the heavy 
nuclei image of Virgo, deflected in the GMF. If one also 
adds the condition that the reconstructed source should be 
located at less than 10° from M87, the above probability 
falls to 2± 3 X 10^^. One may argue that this detection 
can have been triggered by the overdensity of events in the 
"Cen A region" which may be due to another reason. In 
such a case, the relevant probability is the probability to re- 
construct the source at less than 10° from M87 in a data 
set already containing a feature as significant as the Auger 
feature. It is equal to ~ 4 x 10~^. 

4 Discussion 

We first cross-check the result of Section 3 with a "blind- 
like" analysis by dividing the Auger data set in two: The 
first 27 events and the 69 — 27 = 42 newer ones. For the 
first set, we determine the "best sector" on the sky in which 
the significance of the signal is the largest [1], compared 
to some random background made of 27 events following 
Auger exposure and spectrum. This sector is highlighted 
in orange in Fig. 3. It contains 10 out of the 27 events - 
see upper panel. Then, the newer 69 — 27 = 42 events are 
analysed with this sector. It contains 5 events out of 42, and 
Con:(X',l/E) ~ 0.38 -see lower panel. The probability 
to have in some background made of 42 events at least 5 
events in this sector and Corr(X',l/i?) > 0.38 is -- 2%. 
This value is compatible with the order of magnitude of the 
probability computed in the previous Section (0.4%). 

We suggest in Section 3 that the events in the Cen A region 
may be a hint at the heavy nuclei image of Virgo. Let us 
note that the anisotropy above 55 EeV in the Auger data is 
due to the overdensity of events in the Cen A region. By 
computing the 2, 3 and 4-point autocorrelation functions, 
defined as in Ref. [2], we find that the rest of the sky is 
currently still compatible with isotropy. We plot in Fig. 4 
the probability that the signal in the data is a fluctuation of 
the background as a function of the angle (labels "69"). The 
3 and 4-point autocorrelation functions have a minimum on 
the ^ 20° scale. If one removes the events in the 20° region 
around Cen A ("no Cen A"), the minimum disappears. 

Let us also point out that some GMF configurations are 
compatible with the possibility presented in Section 3. 
Fig. 5 shows the iron image, for E = QO — 140 EeV, of 
the Virgo cluster in a reshaped version of existing GMF 
models, whose parameters are discussed in [1]. If one adds 
a turbulent component to this model, low energy events 
would be spread in the whole Cen A region. If Virgo would 
be confirmed in the future to be a heavy nuclei source, it 
would put strong constraints on the GMF. For instance, the 
142 EeV event would give an estimate of UHECR deflec- 
tions in the Galactic Northern halo, and the contribution of 
the turbulent GMF to deflections should be small because 
of the small spread of the image in the Cen A region. It 
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Figure 3: Blind-like analysis with the two consecutive 
Auger data sets (27 and 69 — 27 = 42 events). "Best sec- 
tor" in orange, and red cross for its vertex. Same key as in 
Fig. 2. Upper panel: First data set. The sector contains the 
10 events surrounded with magenta circles; Lower panel: 
Second data set. The sector contains 5 events. 
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Figure 4: Probability that the signal in the 2, 3 and 4-point 
autocorrelation functions is a fluctuation of the background 
as a function of the angle S, both for the sets of 69 events 
(labels "69") and of 56 events, excluding the 20° region 
around Cen A (labels "no Cen A"). 

would also bring constraints on the GMF strengths, exten- 
sions and geometries in the Galactic disk and halo [ 1 ] . 

This would put restrictions on the source(s), too. Due to 
their equal rigidities E/Z, 60 — 80 EeV iron nuclei are de- 
flected as 2 — 3 EeV protons in cosmic magnetic fields. 
Thus, if UHECR sources accelerate both nuclei and pro- 
tons and if the Cen A region events are heavy nuclei, one 
should have protons at Z times lower energies, in the same 
region [21]. Since Auger data is compatible with isotropy 
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Figure 5: Image of UHE iron nuclei emitted by Virgo 
(black disk) and deflected in the regular GMF model dis- 
cussed in [1]. Shaded areas represent the arrival directions 
at Earth of cosmic rays with given energies, see key. 



at low energies, this would imply either that the source 
spectrum is harder than a E^"^ spectrum or that the ratio of 
accelerated protons to nuclei in the source(s) is not larger 
than one to one. The source(s) may for instance have a hard 
spectrum as in the model of Ref. [22]. This would also put 
restrictions on the maximum acceleration energy, as well 
as on the physical conditions in the Virgo cluster [1]. 

Only next generation experiments, such as JEM-EUSO, 
will have enough statistics to confirm or rule out this pos- 
sibility. Auger will triple its statistics during its lifetime. If 
Virgo is the source of the 142 EeV event, it may detect an- 
other event above 10^*^ eV in the same region, which would 
be a hint. It will also settle if the Cen A overdensity is 
not a statistical fluctuation. If it is not a fluctuation, and 
if the Auger composition is correct, there are three main 
possibilities: this overdensity may be due to magnetic lens- 
ing [9, 11], to nuclei from Virgo as proposed above, or to 
nuclei emitted by Cen A [23]. However, Cen A may not be 
powerful enough to be the source [24] . 

5 Conclusions and perspectives 

We have proposed here both a new method to detect some 
UHE heavy nuclei sources in the data, and a new interpre- 
tation of the Auger data which is both consistent with its 
composition and anisotropy measurements. 

In Section 2, we showed that some heavy nuclei sources 
located in some parts of the sky may have at least one of 
their images which looks like an enlarged proton source 
image. We presented a method which is able to detect 
such favourable sources on top of background. In Sec- 
tion 3, we applied this method to the 69 Auger events with 
£: > 55 EeV of Ref. [12] and found that the 142 EeV event 
and the Cen A region events may be the UHE heavy nuclei 
image of the Virgo cluster The associated reconstructed 
source is indeed located at only a few degrees from Virgo, 
at ~ 8.5° from M87. The probability to reconstruct the 
source at less than 10° from M87 in some data already 
containing a comparable set of events is ~ 0.4%. A better 
knowledge of the GMF than currently available, and more 



UHECR data are however still needed to confirm or rule 
out this possibility. 

In Section 4, we cross-checked this possibility with a 
"blind-like" analysis. We also pointed out that the 
anisotropy in the Auger data is due to the Cen A region 
events. We note that the next generation of UHECR exper- 
iments will be able to test the possibility presented here. 
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